Metal-insulator transition in a generalized 
Hubbard model with correlated hopping at 

half-filling 

L. Didukh and V. Hankevych 

Ternopil State Technical University, Department of Physics, 
56 Rus'ka Str., Ternopil UA-282001, Ukraine; Tel. .-+380352251946, 
Fax: +380352254983; E-mail: didukh@tu.edu.te.ua 

Abstract 

In the present paper metal-insulator transition is studied in a generalized 
Hubbard model with correlated hopping at half-filling and zero temperature. 
Single-particle Green function and energy spectrum of electron system are 
calculated. The expressions for energy gap width and the concentration of 
polar states (holes or doublons) are obtained. The conditions for metallic and 
insulating states are found. 

PACS number(s): 71.30.+h, 71.10.Fd, 71.27.+a, 71.28.+d, 

1 Introduction 

In the present paper we consider a generalized Hubbard model including correlated 
hopping |], @j. In such a model the hopping integrals, which describe hoppings of 
holes and doublons are different. These hopping integrals also are distinguished 
from the hopping integral which is connected with the processes of paired creation 
and destruction of holes and doublons. In recent years the similar models have been 
studied intensively [[J - [jnj. In particular, some of these models || [|, [7], have 
been solved exactly under the condition that the number of doubly occupied sites 
is conserved. 

The important puzzle arising in an investigation of the generalized models is 
metal-insulator transition problem. In papers |2j, [L6[ a new mean-field approximation 



(MFA) which leads to the correct description of metal-insulator transition (MIT) has 
been proposed. In the present paper we use this MFA to study MIT in a generalized 
Hubbard model with correlated hopping at half-filling and zero temperature. In 
Sec. 2 we introduce the Hamiltonian of narrow-band model. The decoupling scheme 
of Green functions is described in Sec. 3. Also single-particle Green function and 
energy spectrum of electon system are calculated. In Sec. 4 the expression for energy 
gap width is found. With the help of this formula and the obtained expression for 
the concentration of nolar states MTT is studied. Finallv. Sec. 5 is devoted to the 



2 Hamiltonian of narrow-band model 



Theoretical analysis, on the one hand, and available experimental data, on the other 
hand, point out the necessity of the Hubbard model |l7j generalization by taking 



into account interelectron interactions describing intersite hoppings of electrons (cor- 
related hopping) m, [18], [19| . The characteristic property of these correlated hopping 
integrals is the dependence on the occupation of sites by electrons. 

So, we start from the following generalization of the Hubbard model including 
correlated hopping [JTJ, |2|: 

H = -H afa a iv + H! a t ( *0 + H J ( ik 3 k ) n k J dja 

ia ija \ k / 

+[/£> iT n ii , (2.1) 

i 

where af al ai a are the operators of creation and destruction for the electrons with 
spin o (p =|, |) on z-site, n ia = af a a ia , = + n^; /i is the chemical potential; 
U is the intra-atomic Coulomb repulsion; i is the matrix element which describe 
the hoppings of electrons between nearest-neighbor sites of lattice in consequence of 
electron-ion interaction, 

2 

J{ikjk) = j J 0*(r - R,)0(r - R^^L— |0( r ' - R fc )| 2 drdr', (2.2) 

(0-function is the Wannier function). The prime at second sum in Eq. (2.1) signifies 
that i ^ j. 

In Hamiltonian (2.1) we rewrite the sum Y^jak J{^3 k ) a tr n k a ja in the form 
Y J {ikjk)af a n k a j(7 + ^ (j (iiij)af a a jrT n^ + h.c.^j (2.3) 

ija k^jti ija 
k^j 

(a denotes the spin projection which is opposite to a), here we have used that 
J(iiji) = J(jiii) = J(iiij) in consequence of the matrix elements symmetry. Let us 
suppose (as in the papers []IJ, Q) that 

J2J2 J ( ik J k ) a ta n kaja = T^'at a ja (2.4) 

ija k^i ija 
k+j 

with T% = nj^k^i J{ikjk) and n = (n^ + na) (sites i and j are nearest neighbors); 

k+j 

it should be noted that this supposition is exact in the homeopolar limit (n» = 1). 
Thus at half- filling (n = 1) we can write Hamiltonian (2.1) in the form 



H = -/j, Y <4a a iv + t Y^ a ta a jcr + (af a a ja niv + h.c.^ 

ia ija ija 

+Uj2^nn, (2.5) 



where t = t + T\ and T 2 = J(iiij). 



dependence of the hopping integral t. It is the distinction of the present model from 
similar models || |], ||, |7], |18|, |19| . In the case n = 1 taking into account the correlated 



hopping T\ leads to the renormalization of the hopping integral t only. 

Rewrite Hamiltonian (2.5) in terms of X^-Hubbard operators |2(| using the 
formulae 

a i1 — A i — A i j a *T — A i — A i j 

a il ~ A « + A » ' a il ~ A i + A i ' 

where X^ z is the transition-operator of i-site from state \l) to state |/c); |0) denotes 
the site, which is not occupied by an electron (hole), \a) = af a \0) denotes the singly 
occupied (by an electron with spin cr) z-site, |2) = a^a^O) denotes the doubly 
occupied (by two electrons with the opposite spins) i-site (doublon). 
In terms of X^- operators Hamiltonian (2.5) takes the following form: 

H = H + H 1 + H[, (2.6) 

with 



Ho = -fiJ2 (4 + *i + zx? ) + uj^xi 

i i 

H x = t^xfxf + igxfxf, 

ija ija 

H[ = t'J2' (r)«X°°Xp + h.c.) , 



ija 



where Xf = Xf l X[ k is the operator of the number of |/c)-states on z-site, rj^ = 

-i, vi = l; 

t = t + 2T 2 , t' = t + T 2 . (2.7) 
The single-particle Green function 

G^(E) = «a p >+\) (2.8) 
in terms of Hubbard operators is written 

g;A e ) = ((x° P a \x;P}) + vA(x°/\x 2 /)) + n a ((xp |x;,°» 

+ ((X?\X 2 /}}. (2.9) 
Hq describes the atomic limit of narrow-band models. 

Hi describes the translational hopping of holes and doublons. In the present 
model (in contrast to the narrow-band models of the Hubbard type) the hopping in- 
tegrals of holes t and doublons t are different. It should be noted that in consequence 
of the difference of the hopping integrals, which describe translational hopping of 
current carriers within the lower (hole) band and upper (doublon) band, the energy 
width of the upper band can be much smaller and the effective mass of current 
carriers within this band can be much larger than in the lower band. Thus, within 
the proposed model the ideas of the "wide" and "narrow" subbands and the "light" 
and "heavy" current carriers are introduced (as the result of electron-electron inter- 
actions). 



3 Energy spectrum of electron system: MFA 

Here we use MFA proposed in papers 0, [16] to calculate the energy spectrum of 
electron system described by Hamiltonian (2.7). 

The Green function ((Xp a \XpP)) is given by the equation 



{e + p)((x?\xp>)) = 6 -^(x; + x° p ) + (([x /,^] |x;,°» 



+{{x»/,h{ \x?)) 



with [A, B] = AB — BA, 



yOu tt 



t J2 ((x; + x Q p )xf + x™xf ) - tY^xfxf, 



(3.1) 



(3.2) 



v0& TJl 



era \ra2 
j 



-t']Tx p 02 A7 + t'Y J x; a x 



(3.3) 



To break off the sequence of Green function equations according to generalized 
Hartree-Fock approximation we suppose that 



yOff tt 
A-p , -Hi 



£e(R7)*f, X^,H[ = $>(pj)Aj 



(3.4) 



where e(pj) and e\{pj) are the non-operator expressions. The representation choice 
of the commutators in form (3.2) and (3.3) is prompted by the operator structure of 
these commutators which maps the energy non-equivalence of the hopping processes 
prescribed by Hi and H[. Taking into account (3.4) we rewrite Eq. (3.1) in the form 

(e + ^)((x°/\x;P)) = 5 -f{x; + xD + E^')((^I^; )) 

3 

+ E^to)((^f|x;, )). (3.5) 



After anticommutation of both sides of the first of formulae (3.4) with X% and the 



second formula with X% a we obtain 



e{pk){Xl + X i 



GO 

V 



t{x; + x Q p ){x% + x° k ) + txfx i 

Spkt £ X% Xj a + 5 pk i^2 XfX^ 



Zv20 y02 

-tA k A p , 



(3.6) 



e 1 (pk)(X° +Xl 



-t\x; + X${Xl + Xl) + t'X™X] 
-5 p kt' ^2 xj°x® a + 5 pk t' xi a xj 



Similarly, for the Green function ((Xp 2 \XpP)) we can write the equation 

(e + ^-u)((x; 2 \x;P)) = J2~e(pj)((xp\x° }) 

j 

+E e 2 to)((^n^ )), (3.8) 

j 

where e(pj) and ^{pj) are determined through the expressions which are analogous 
to (3.6) and (3.7). Thus we obtain the closed system of equations for the Green 
functions ((X^\X°P)) and ((X™ \X;P)). 

By neglecting correlated hopping and by averaging expressions (3.6) and (3.7) 
we obtain the approximations |17|, ^3], [24|; the defects of these approximations are 
well-known (see, for example Ref. P5| ). Here we use the approach which has been 
proposed in the papers 0, pL6| . 

To determine e(pj), ei(pj) we rewrite Xf'-operator in Eqs. (3.6) and (3.7) in the 
form |26| Xf l = af k au, where af k) an are the operators of creation and destruction 
for | A;)- and |/)-states on i-site respectively (the Schubin-Wonsowsky operators [pTfl); 
thus X® = a^a i0 , X 2 = a^a^, X[ = af a a ia . Let us substitute a-operators by c- 
numbers in Eqs. (3.6) and (3.7) (here there is a partial equivalence with slave boson 



method 28 



OLicj = , a i0 = a i0 = a i2 — a i2 



d 1 ' 2 (3.9) 



(we consider a paramagnetic case, electron concentration on site n = 1); d is the 
concentration of polar states (holes or doublons). 

The proposed approximation is based on the following physical idea. Let us 
consider a paramagnetic Mott-Hubbard insulator at T ^ 0. Within the wide tem- 
perature interval (fc^T -C U) the concentration of polar states is small (d <C 1). An 
analogous consideration is valid for a paramagnetic Mott-Hubbard semimetal (hole 
and doublon subbands overlap weakly, (i < 1). So, the change of states and polar 
excitations influences on |<r)-states weakly. Thus we may consider |<r)-states as the 
quasiclassical system and substitute the operators af a , a ia by c-numbers. In addi- 
tion, when we find e(pj), ei(pj) we substitute the creation and destruction operators 
of |0)- and |2)-states through the respective quasiclassical expressions. Actually the 
proposed approximation is equivalent to a separation of the charge and spin degrees 
of freedom. Note that the present approach is the most justifiable when d —>■ 0. 

Thus in k-representation we obtain 

e(k) = (1 - 2d + 2d 2 )t k - 2d%, Cl (k) = -2dt' k , (3.10) 

where t^, 4, are the Fourier transforms of the hopping integral t, t, t' respec- 
tively. Similarly, we find that 

e(k) = (1 - 2d + 2d 2 )t k - 2d%, e 2 (k) = -2dt' k . (3.11) 

The Fourier transform of the Green function ((Xp a \XpP)) is found from the 
system of equations (3.5) and (3.8) 



with 



ElAk) = ., + (l^tM±E^ Fk , (3 ,3) 

F k = sj [S(i k - f k ) - [/] 2 + (4<ft' k ) 2 , 5 = 1 - 2d + 4d 2 . (3.14) 

An analogous procedure is realized also in the equations for the other Green 
functions in Eq. (2.6). 

Finally, in k-representation the single-particle Green function is 

^ = s(^oo + F^«)- (3 ' 15) 

A^-- 2 f , * = i + ^. (3.16) 

Single-particle Green function (3.15) gives the exact atomic and band limits: if 
U — and i k = t k = t k = to(k) (it means neglecting correlated hopping) then 
G ^(E) takes the band form (d — 1/4 when U = 0), if i k = t k = t k — > then we 
obtain the exact atomic limit. 

The peculiarities of obtained quasiparticle energy spectrum (3.13) of narrow- 
band system which is described by Hamiltonian (2.5) are the dependence on the 
concetration of polar states and the non-equivalence of the lower and upper Hubbard 
bands. This non-equivalence is caused by the difference of the hopping integrals t, 
t, t'. 

Quasiparticle energy spectrum (3.13) allows to study MIT in the proposed model. 



4 Metal-insulator transition 

With the help of energy spectrum of electrons (3.13) we find the expression for the 
energy gap width (difference of energies between bottom of the upper and top of 
the lower Hubbard bands): 

AE = -(1 - 2d)(w + w) + h < Q 1 + Q 2 ), (4.1) 

Qi = ^[B(w-w)-Uf + (4dzt>) 2 , 
Q 2 = ^[B(w-w) + U] 2 + (4dzt') 2 , 

where w and w are the half-widths of the lower (hole) and upper (doublon) Hubbard 
bands respictevely: w = z\t\, w = z\t\ (z is the number of nearest neighbors to a 
site) . 

The peculiarities of the expression for energy gap (4.1) are dependences on the 
concentration of polar states, on the widths of hole and doublon bands, on the hop- 
ping integral t' (thus on external pressure). At given U, t, i, t' (constant external 
pressure) the concentration dependence of AE allows to study MIT under the ac- 
tion of external influences: temperature change, photoeffect and magnetic field. In 



transition from the state of a paramagnetic metal to the paramagnetic insulator 
state in the (Vi-zCr^O^ compound 0, [H], in N1S2 |[2) and in the N^-iSe^ 
system |33] should be noted). Under the action of light or magnetic field the 

concentration of polar states can be changed; it leads to the fact that the energy 
gap width is changed also and MIT can occur. 

Distinction of formulae (3.13)-(3.16), (4.1) from earlier obtained results (e.g., 
see reviews 



31| , |3q , |36fl ) is the dependence on concentration of polar states. Let us 



find the expression for its calculation. 

The concentration of polar states is given by the equation 



d = (X?) 



+00 

/ ME)dE 




exp^M + 1 



+ 



exp^M + 1 



(4.2) 



where 



B(t k - t k ) U 



2F k 



2F k 



1 ggkj^k) u 

2 2F k 2F k : 



9 = kgT, kg is the Boltzmann's constant, N is the number of sites, J k (E) is the 
spectral intensity of the Green function 



1 



«*W»k = T~ 



c k 



+ 



D k 



4tt {E-E^k) E-E 2 (k) 
At T = and the rectangular density of states 



(4.3) 



— y s(E-t k ) = —i 

k 



(w 2 - E 2 



(9(x) — 1 if x > 0, =0 otherwise) from Eq. (4.2) we obtain that 



B t — t U 
[(p(e ) - <^(-e )J + 



B 2 (i-tf 



x In 



z X 



Vx 



x 



Ae - BU{t - t) 



y^X(p(-e ) + Ae - BU(t - t) 



with 



fill — fi 2 



U-2 (U <w + w) 



'1 - 2d + 2d 2 )w - 2d 2 w 



(4.4) 



'\| (l-2d) 2 (t + t) 2 -A' r (l-2d)(w + w 

(p(e) = {\e 2 -2BU(t-t)e + U 2 y\ X = B 2 (t — t) 2 + (Adt') 2 
For narrow-band semimetal (d <C 1) Eq. (4.4) takes the following form: 



U. 



Fig. 1 shows the dependence of d on U/w which is obtained from Eq. (4.4). The 
parameters t\ = Ti/|t |, t 2 = T 2 /\to\ characterize the value of correlated hopping. 
One can see that a value of d depends on the parameters of correlated hopping t% , r 2 
(thus on w/w) weakly when U/w is close to zero. But with the increase of U/w the 
concentration of polar states becomes strongly dependent on the parameters T\, r 2 . 
It testifies on the fact that taking into account the correlated hopping is important 
to consider the metal-insulator transition problem. 

Fig. 1 shows also that if U > w + w then the concentration of polar states d = 0. 
In the special case t + 1 = t' = this consequence is in accordance with the results 
of Refs. |, 0,0. 

At T = the energy gap width AE < (i.e. MIT occurs) when the condition 
U<w + w (4.6) 



is satisfied (in agreement with general physical ideas fl3l| ). For the special case t' = 
condition (4.6) covers the exact results of Refs. [|, |, §]. 

Fig. 2 which is obtained from formula (4.1) using Eq. (4.4) shows that in a 
metallic state the overlapping of energy subbands decreases and in an insulating 
state the energy gap width increases with decrease of the parameter w/w (at given 
U/w). 

In the Hubbard model energy gap width (4.1) takes the following form: 



AE = -2w(l - 2d) + ^U 2 + (4dw) 2 } (4.7) 
and the concentration of polar states (4.4) is 

d =Q + 3^ in ( i - 4d )) ,,(2t "- c,) - (4 - 8) 

In the region of metal-insulator transition d = 1/4 — U/(8w); this dependence is in 
qualitative accordance with the result of Brinkman and Rice [37] obtained by use of 



Gutzwiller variational method [38|, those of the general Gutzwiller-correlated wave 



functions in infinite dimensions and the Kotliar-Ruckenstein slave bosons 
For U/2w — > we obtain d = 1/4 + U/(8w) ln(U/2w) (if we consider Coulomb 
repulsion as perturbation then d(U — > 0) = 1/4 — 0(U)); in order to compare the 
obtained dependence (4.8) d on U/w in the Hubbard model with other approximate 
theories see e.g. [0). AE < when the condition 2w > U is satisfied. 



5 Conclusions 

In the present paper Mott-Hubbard transition has been studied in a generalized 
Hubbard model with correlated hopping at half-filling using a mean-field approxi- 
mation ||, |i~6fl . 

We have obtained the expression to calculate the concentration of polar states. 
With the increase of the parameter w/w (0 < w/w < 1) the concentration of polar 
states increases at given U/w. 

Quasiparticle energy spectrum has been calculated. With the help of this en- 



non-equivalence of the upper and lower Hubbard bands. The increase of the param- 
eter w/w (at given U/w) leads to decreasing energy gap width. In consequence of 
it, in particular, MIT occurs when the condition U/w — 1 + w/w is satisfied. 

The cases n ^ 1 and T ^ 0, an application of the obtained results to the 
interpretation of the experimental data will be considered in the next papers. 

The authors are grateful to Prof. I. V. Stasyuk for valuable discussions. 

References 

[1] L. Didukh, Fiz. Tverd. Tela 19, 1217 (1977); L. Didukh, Ukrainian-French 
Simposium "Condensed Matter: Science and Industry". Abstracts, p. 275 (Lviv, 
1993). 

[2] L. Didukh, Journ. of Phys. Stud. 1, 241 (1997). 

[3] R. Strack, D. Vollhardt, Phys. Rev. Lett. 70, 2637 (1993). 

[4] A. A. Ovchinnikov, J. Phys.: Condens. Matter, 6, 11057 (1994). 

[5] J. E. Hirsh, Physica B 199-200, 366 (1994). 

[6] A. Schadschneider, Phys. Rev. B 51, 10386 (1995). 

[7] J. de Boer, V. E. Korepin, and A. Schadschneider, Phys. Rev. Lett. 74, 789 
(1995). 

[8] A. A. Aligia, L. Arrachea, and E. R. Gagliano, Phys. Rev. B 51, 13774 (1995). 

[9] M. Airoldi and A. Parola, Phys. Rev. B 51, 16327 (1995). 

[10] H. Q. Lin and J. E. Hirsch, Phys. Rev. B 52, 16155 (1995). 

[11] E. R. Gagliano, A. A. Aligia, and L. Arrachea, M. Avignon, Phys. Rev. B 51, 
14012 (1995). 

[12] E. R. Gagliano, A. A. Aligia, L. Arrachea, and M. Avignon, Physica B 223-224, 
605 (1996). 

[13] G. Santoro, N. Manini, A. Parola, and E. Tossatti, Phys. Rev. B 53, 828 (1996). 
[14] L. Arrachea, E. R. Gagliano, and A. A. Aligia, Phys. Rev. B 55, 1173 (1997). 
[15] B. Chattopadhyay and D. M. Gaitonde, Phys. Rev. B 55, 15364 (1997). 
[16] L. Didukh, Phys. Stat. Sol. (B) 206, R5 (1998). 

[17] J. Hubbard, Proc. Roy. Soc. A 276, 238 (1963); J. Hubbard, Proc. Roy. Soc. 
A 281, 401 (1964). 

[18] S. Kivelson, W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. Lett. 58, 
1899 (1987); 60, 72 (1988); D. Baeriswyl, P. Horsch, and K. Maki, Phys. Rev. 



[19] J. E. Hirsch, Physica C 158, 326 (1989); J. E. Hirsch and F. Marsilio, Phys. 
Rev. B 41, 2049 (1990); F. Marsilio and J. E. Hirsch, Phys. Rev. B 41, 6435 
(1990). 

[20] J. Hubbard, Proc. Roy. Soc. A 285, 542 (1965). 

[21] L. D. Didukh, I. V. Stasyuk, Ukr. Phys. Journ. 13, 1262 (1969). 

[22] L. M. Roth, Phys. Rev. Lett. 20, 1431 (1968); D. N. Zubarev, Yu. G. Rudoj, 
Uspekhi fizicheskikh nauk 163, 103 (1993). 

[23] A. B. Harris and R. V. Lange, Phys. Rev. 157, 295 (1967). 

[24] L. Roth, Phys. Rev. 184, 451 (1969). 

[25] A. Kawabata, Progr. Theor. Phys. 48, 1793 (1972). 

[26] L. Didukh, I. Stasyuk, Fiz. Met. Metalov. 26, 435 (1968). 

[27] S.Schubin, S.Wonsowsky, Proc. Roy. Soc. A 145, 159 (1934). 

[28] G. Kotliar and A. E. Ruckenstein, Phys. Rev. Lett. 57, 1362 (1986). 

[29] L. Didukh, V. Hankevych, in: INTAS-Ukraine Workshop on Cond. Mat. Phys., 
Abstracts, p. 105, Lviv (1998). 

[30] D. B. Mc Whan, J. B. Remeika, Phys. Rev. B 2, 3734 (1970). 

[31] N. F. Mott, Metal-Insulator Transitions Taylor and Francis, London 1990. 

[32] J. A. Wilson, The Metallic and Nonmetallic States of Matter, Taylor & Francis, 
London 1985. 

[33] X. Yao and J. M. Honig, T. Hogan and C. Kannewurf, J. Spalek, Phys. Rev. B 
54, 17469 (1996). 

[34] X. Yao, Y.-K. Kuo, D. K. Powell, and J. W. Brill, J. M. Honig, Phys. Rev. B 
56, 7129 (1997). 

[35] Yu. A. Izyumov, Uspekhi Fizicheskikh Nauk 165, 403 (1995). 

[36] F. Gebhard, The Mott Metal-Insulator Transition - Models and Methods, 
Springer, Berlin 1997. 

[37] W. F. Brinkman and T. M. Rice, Phys. Rev. B 2, 4302 (1970). 

[38] M. C. Gutzwiller, Phys. Rev. Lett 10, 159 (1963). 

[39] F. Gebhard, Phys. Rev. B 41, 9452 (1990). 

[40] M. Dzierzawa, D. Baeriswyl, and M. Distasio, Phys. Rev. B 51, 1993 (1995). 

Figure 1: Concentration of polar states d as a function of U/w: the upper curve 
corresponds to T\ = t 2 = 0; the middle curve - T\ = r 2 = 0.1; the lower curve - 

n = T 2 = 0.2. 



